Novel neutron resonance mode in (1^2 _y2 superconductors 
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We show that a new resonant magnetic excitation at incommensurate momenta, observed recently 
by inelastic neutron scattering experiments on YBa2Cu3 06.85 and YBa2Cu306.6, is a spin exciton. 
Its location in the magnetic Brillouin zone and its frequency are determined by the momentum 
dependence of the particle-hole continuum. We identify several features that distinguish this novel 
mode from the previous resonance mode observed near Q = (7r,7r), such as its intensity maximum 
which occurs in a different part of the magnetic Brillouin zone. 

PACS numbers: 71.10.Ca,74.20.Fg,74.25.Ha,74.72.Bk 



While it seems established that the resonance peak is 
a universal feature of the high-temperature superconduc- 
tors 0,0,0, its origin, its role for the pairing process, 
and the effects arising from its interactions with elec- 
trons are still intensively debated 1^. The peak's inten- 
sity is the highest at Q = (tt, tt) and in YBa2Cu306+x 
(YBCO), where it was studied in great detail, its fre- 
quency QresiQ) follows the Same doping dependence as 
Tc, with QresiQ) ^ 41 meV near optimal doping. As 
one moves away from Q, the peak disperses downwards 
and its intensity decreases rapidly, vanishing around 
Qo = (O.Stt, O.Stt). The doping dependence of Qres{Q)^ 
the downward dispersion of the resonance, and the fact 
that Qo coincides with the distance between nodal (diag- 
onal) points on the Fermi surface are consistent with the 
theoretical idea that the resonance peak is a particle-hole 
bound state below the spin gap (a spin exciton)A] (for 
a review of other theoretical scenarios, see Ref. p). 

Recent inelastic neutron scattering (INS) experiments 
on YBa2Cu306.85 7] superconducting (SC) state 

detected a new resonant magnetic excitation at in- 
commensurate momenta, but at frequencies larger than 
^res(Q)- A similar result was obtained for underdoped 
YBa2Cu306.6 .8]- This new resonance mode is partic- 
ularly pronounced along the diagonal of the magnetic 
Brillouin zone (MBZ) at q ^ Qo- It was suggested 
that this new resonance is a particle-hole bound state 
with an upward dispersion originating at Q (see Fig. 5a 
in Ref. 0). 

In this Letter, we show both analytically and numer- 
ically that the new resonance mode is indeed a spin ex- 
citon that emerges below Tc due to a feedback effect on 
the collective spin excitations arising from the opening 
of the SC gap in a (1^2 _y2 superconductor. We demon- 
strate that the new resonance appears only at momenta 
less than Qo, and is separated from the previous reso- 
nance by a region near Qo in which no resonance exists 
(the "silent band" of Ref. jlj). Thus, the new resonance 
does not form an upward dispersing branch originating 



at Q. We identify several qualitative features that dis- 
tinguish this new resonance (the Q* mode) from the old 
one (the Q mode). In particular, we show that while the 
intensity of the Q mode is largest along q = (tt, t^tt) and 
q = {rjTT^Tr), the Q* mode has its largest intensity along 
the MBZ diagonal. 

We begin by presenting our numerical analysis of the 
new Q* mode. Its emergence can be understood within 
an RPA approach for which the spin susceptibility is 
given by 



x(q,^) = 



xo(q,^) 



1 -^(q)xo(q,^) 



(1) 



where ^(q) is the fermionic four-point vertex, and 
Xo(q7^) is the free-fermion susceptibility, which in the 
SC state is given by the sum of two single bubble dia- 
grams consisting of either normal or anomalous Greens 
functions Jsl, 1^. For our numerical calculation of 
Xo(q, ^), we used a SC gap with dx2_y2 symmetry and a 
normal state tight binding dispersion 

ek = — 2t (cos kx + cos ky) — At' cos kx cos ky — ji (2) 

with t = 250 meV, t' /t = -0.4, and ii/t = -1.083. The 
Fermi surface (FS) obtained from Eq. Q is shown in 
Fig. |3Ja). It describes well the FS measured by photoe- 
mission experiments on Bi2Sr2CaCu208+(5 jlC^ . 

Our main results are presented in Figs. and [21 in 
which we plot Imx(q, oo)^ obtained from a numerical eval- 
uation of Eq.Q, along q = /^(tt, tt) in the SC state. The 
intensity plot of Imx(q, Cc;) shown in Fig. possesses all 
salient features observed in the INS experiments. First, 
we identify a downward dispersion of the Q mode (indi- 
cated by a red arrow). Second, the Q* mode (indicated 
by a yellow arrow) is located at frequencies larger than 
the frequency of the Q mode at (7r,7r), and is confined 
to a small region of momentum space near Qo. The mo- 
mentum position of the Q* mode is almost independent 
of energy. Third, there exists a region in momentum 
space around Qo that separates the Q mode from the Q* 



2 




FIG. 1: (color) RPA results for magnetic excitations in a 
dx2_y2 superconductor. Imx obtained from Eq. (QJ as a func- 
tion of momentum (along q = ?7(7r,7r)) and frequency in the 
SC state. We used Ak = Ao(cos/ccc — cos ky)/ 2 with Ao = 42 
meV, and ^(q) = go[l — 0.1(cos^cc + cos^^)], with go = 0.573 
eV in order to reproduce the correct energy position of the 
Q* mode near 0.8(7r,7r) and the Q mode at (tv^tv). 

one (the "silent band" of Ref. 't^). Note, that as Qo is 
approached from (7r,7r), the Q mode frequency, Qres{^)^ 
as well as its intensity, rapidly decreases. In Fig. [2 we 
present Imx along q = /^(tt, tt) for several frequencies. We 
clearly see that the two modes are separated in momen- 
tum and frequency space. The shaded area in Fig.[21rep- 
resents the silent band, in which Imx is strongly reduced 
from its resonance values. We also find that the position 
of the Q* mode is almost frequency independent, with a 
maximum intensity at ~ 54 meV, in agreement with 
the experimental observations (see Fig. 2a of Ref. ^]). 
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FIG. 2: Im x(Q5^) along q = r]{7v,7v) for several frequencies. 
The position of the silent band is indicated by the shaded 
area. These results reproduce the experimental INS data, see 
Fig. 2(a) in Ref. Jj. 



In what follows, we discuss the physical origin of our 
results. We begin by reviewing the emergence of the Q 
mode. Its origin lies in a feedback effect on the spin 
excitation spectrum (i.e., on xo) due to the opening of 
the SC gap. This feedback effect is universal in that it 
only requires that the FS possesses points, k and k', that 
can be connected by Q (so-called "hot spots") and that 
sgn(Ak) = — sgn(AkO- To demonstrate the universality 
of the resonance mode, we rewrite the denominator of 
Eq. (QJ as 

eQ'-^(Q)Axo(Q,^^) (3) 

where = I - 9{Q)x^'{Q.O) and Axo(Q,^^) = 

XoCQ, ^)~Xo^{Q^ 0)- The main contribution to the static 
normal state susceptibility Xo^iQ^^) comes from high 
internal frequencies in the fermionic bubble. This im- 
plies that Xo^ cind thus are non-universal quantities 
that depend on the details of the band structure. How- 
ever, the appearance of the resonance mode only requires 
> 0, which is satisfied in the paramagnetic state. In 
contrast, Axo(Q,^) is universal since its main contribu- 
tion comes from small internal frequencies of order A. 
Axo can therefore be evaluated by simply linearizing the 
fermionic dispersion near k and k^ The integration over 
momentum llj yields at T = 

where u;± = a; ± 1^/2, and the summation runs over 
all pairs of FS points k and k' separated by Q. There 
are eight fermionic scattering processes from k to k' in 
the first zone. The momenta k involved are those in 
which the boundary of the magnetic zone (defined by 
coskx = — cos ky) crosses the FS. Two of these scatter- 
ing processes are direct with k — k' = Q, two involve 
umklapp scattering with k — k' = Q — (27r, 27r), and four 
involve umklapp scattering with k — k' = Q — (27r, 0) and 
k — k' = Q — (0, 27r). A pair of FS points that are con- 
nected by direct scattering via Q is shown in Fig. OJa) 
(Q is represented by a dashed arrow). 

In the normal state, all eight processes equally con- 
tribute to Axo(Q,^), and one obtains from Eq. 
Axo(Q,^) = -ilQ^^ which identifies with the 

Landau damping rate In the superconducting state, 
Axo becomes a complex function. According to Eq. Q), 
its imaginary part vanishes below a critical frequency 
^c(Q) = |Ak| + |Ak'|, which is the same for all eight 
scattering channels. The dx2_y2 symmetry of the SC gap 
implies Ak = — Aks resulting in a discontinuous jump 
of Imxo at l^c(Q) from zero to 7r7Q'^|Ak| 0. Simulta- 
neously, ReAxo(Q, ^) > is non-zero, diverges logarith- 
mically at Vtc{Q)i and scales as Q? at small frequencies. 



3 



ReAxo therefore varies between at ^2 = and oo at 
Q = Qc' Since > 0, one finds that for any posi- 
tive ^Q, xiQ^^) [Eq. O] acquires a pole at a frequency 
Qres < ^c(Q) where ReAxo(Q, ^res) = S>q^ / 9q and 
ImAxo (Q, ^res) = 0. Imx thus exhibits a ^—function 
at r^res, representing a spin exciton below the particle- 
hole (p-h) continuum. 




FIG. 3: (a) FS of Eq. Q and magnetic scattering vectors, (b) 
Momentum dependence of the critical frequencies Qi^^ (i — 
1,2,3) and the direct gap along q = r]{7T,7T) (see text). 
The open squares (circles) represent the position of the Q 
(Q*) resonance. 

We next consider the dispersion of the exciton. For 
q 7^ Q, the degeneracy of the scattering channels is 
lifted. In particular, for momenta along q = 77(7r,7r), 
one now has three different critical frequencies, rt^\q) 
{i = 1,2,3) 5]. ft^^ is the critical frequency associated 

(2) 

with direct scattering, Qc with umklapp scattering in- 
volving q— (27r, 0) and q— (0, 27r), and l^c^^ with umklapp 
scattering by q — (27r, 27r). The scattering momenta for 
q close to Qo are shown in Fig. I^a), and Q^\q) for 
q = 77(7r,7r) are presented in Fig. [3fb). For all of these 
scattering processes, Ak and Ak' still have opposite signs 
jl3j . As a result, ImAxo exhibits three discontinuous 
jumps at Q,^c\^)^ and ReAxo diverges logarithmically 
at Qc\^) 0. However, ImAxo is zero only below the 
smallest (q) , and hence a true resonance is only pos- 
sible below the smallest critical frequency. The splitting 
of the critical frequencies for q 7^ Q was also discussed 
in Ref. Q. 

It follows from Fig. O^b) that for O.S < r] < 1.2, the 
smallest critical frequency, corresponds to direct 

(i.e., non-umklapp) scattering. decreases away from 
T] = 1 and eventually vanishes at r] = 0.8, when direct 
scattering occurs between nodal points at the Fermi sur- 
face 5] . Since the exciton is necessary located below , 
its frequency also decreases and eventually vanishes at 
T] = 0.8. Moreover, upon approaching rj = O.S the jump 
in ImAxo(q, ^c^^) decreases as k and approach the 
nodal points. Accordingly, the resonance frequency ftres 
moves closer to and the intensity of the resonance 

decreases j5|J. 



At Qo = 0.8(7r, tt), Imxo(Qo, ^) is non-zero for (7 > 0, 
and one finds Imx(Qo,^) = 7q^^ at small frequencies, 
where 7q^ = Tq*^!!^ and va is the gap velocity at 
the nodal points. The factor 1/8 arises since only a sin- 
gle (direct) scattering channel contributes to 7q^, while 
eight channels contribute to Jq^ . However, since the 
Fermi velocities at the nodal points are antiparallel, 
depends on vf only through the ratio vf/va ^ 20 which 
compensates the small prefactor 15j. As a result, 7q^ is 
comparable to 70^*^, thus giving rise to a weak and fea- 
tureless frequency dependence of Imx(Qo,^) similar to 
that of Imxo(Q5^) in the normal state. The vanishing 
of the gap in the p-h continuum at Qo together with the 
large value of 7q^ explains the experimental observation 
of a "silent band" in Ref. (the position of which is 
indicated by the shaded area in Fig.EJ. 

For momenta q < Qo, (i.e., 77 < 0.8) the nodal points 
cannot be connected, and a direct gap opens for excita- 
tions into the p-h continuum. This gap is independent of 
the SC gap, and given by = v^? • (Qo — q) (see dashed- 
dotted line in Fig.jSfb)). Due to a large |vi?|, becomes 
equal to Q^"^ close to Qo at q^^ = rj^Q with rjd = 0.773 
(i.e., Qo — qd = 0.027Q for the dispersion of Eq. For 
Vd < V < 0-8 5 finds that Imxo vanishes below 1]^, and 
Imxo — ^ > Hence, Rexo does not di- 

verge at , and no resonance peak exists in this region, 

extending the silent band. However, Rexo possesses a 

(2) 

logarithmic divergence at , and hence it satisfies the 
resonance condition ReAxo(q, ^res) = Cq^/^'q? at some 
frequency fires below Qc • Once Q.^ crosses 
damping at Ctres vanishes and a true pole in Imx oc- 
curs, leading to the appearance of the Q* mode (open 
circles in Fig. Ol^b)). As one moves further away from 
Qo, one finds that the Q* resonance is rapidly suppressed. 
This suppression arises from the rapid decrease of Vtc as 
well as the decrease of the bare static spin susceptibility, 
Xo(q7 0). This behavior of Xq{^.->^) is similar to that in 
a system of free fermions with a circular Fermi surface, 
if we identify Qo with 2pi?. Both effects lead to a rapid 
shift of the Q* mode towards the edge of the p-h contin- 
uum, and to a decrease in its intensity. As a result, this 
resonance is only visible near Qo. 

The Q and Q* modes are not only separated in fre- 
quency, as discussed above, but their intensity maxima 
are also located in different parts of the MBZ; this rep- 
resents a major qualitative distinction between the two 
modes. In Fig. 01 we present intensity plots of Imx as 
a function of momentum for O > Vtres{Q)i probing the 
Q* mode [Fig. 2fb)], and for < ftresiQ)^ probing the 
Q mode [Fig. ^a)]. The difference is striking. While 
the intensity of the Q mode is largest along q = (tt, t^tt) 
and q = (7^7r,7r), the Q* mode has its largest intensity 
along the diagonal direction, i.e., along q = 77(7r,7r) and 
q = [(2 — 77)7r, tjtt]. This rotation of the intensity pattern 
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FIG. 4: (color) Intensity plot of Imx, Eq. (QJ, as a func- 
tion of momentum (blue (red) implies a small (large) value of 
Imx) at: (a) 33 meV (Q mode), and (b) 57 meV (Q* mode). 
The Q mode forms a distorted ring in momentum space, with 
intensity maxima along {"K^q). In contrast, the intensity of 
the Q* mode is largest along the zone diagonal. In (b), the 
periodic structure of Imx around Q reflects the momentum 
dependence of the p-h continuum at this energy. 

by 45° reflects the qualitative difference in the origin of 
the two modes. The intensity of the Q mode is at a maxi- 
mum along q = (tt, t^tt) and q = (t^tt, tt), since in this case 
the fermions that are scattered by q are located farther 
from the nodes than for diagonal scattering. In contrast, 
the Q* mode arises from the rapid opening of a gap in 
the p-h continuum below Qo, which is most pronounced 
along the diagonal directions of the zone. The rotation of 
the intensity pattern by 45°, which is the most significant 
distinction between the Q an d Q * modes, was observed 
in recent INS experiments |3, 

In summary, we showed that a new resonant mag- 
netic excitation at incommensurate momenta, observed 
recently by inelastic neutron scattering experiments on 
YBa2Cu306.85 7] YBa2Cu306.6 .8], is a spin exci- 
ton arising from umklapp scattering. Its location in the 
zone and its frequency are determined by the momen- 
tum dependence of the p-h continuum. It is confined to 
a small region in q just below Qo, and is separated from 
the Q resonance by a "silent band" where Imx is strongly 
suppressed. We also found that the intensity maxima of 
the two modes are rotated by 45° relative to each other. 
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